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CD ; ABSTRACT 

We use a simple toy-model to discuss global MHD modes of a neutron star, taking 
f**«> \ into account the magnetic coupling between the elastic crust and the fluid core. Our 

results suggest that the notion of pure torsional crust modes is not useful for the 
coupled system. All modes excite Alfven waves in the core. However, we also show 
that the modes that are most likely to be excited by a fractured crust, eg. during a 
\Q \ magnetar flare, are such that the crust and the core oscillate in concert. For our simple 

model, the frequencies of these modes are similar to the "pure crustal" frequencies. In 
*-Q \ addition, our model provides a natural explanation for the presence of lower frequency 

(< 30 Hz) quasi-periodic oscillations seen in the December 2004 giant flare of SGR 
1806-20. 



^p. 1 INTRODUCTION 



The recent observational evidence {Israel et aJ kOQfil Istrohmaver fc WattslbOQfih of quasi-periodic oscillations (QPOs) during 
C/j , giant flares in the soft gamma-ra y repeaters ( SGRs ) 1806 — 20 and 1900 + 14 may herald the beginning of an exciting new era 



giant flares in the soft gamma-ra y repeaters ( SGRs ) 1806 — 20 and 1900 + 14 may herald the beginning of an exciting new era 
for neutron star physics (see also lBarat et, all l|l983f ) for older evidence of a QPO in the flare of SGR 0526-66). These extremely 
viole nt events can be understood within the "magnetar" mo del, introduced by Duncan and Thompson more than a decade 
ago iDuncan &: Thompsonlll992l ; I Thomson fc Duncanlll995l) . Magnetars are neutron stars endowed with powerful magnetic 
fields (~ in 14 - 15 Q) masquerading as both SGRs and anomalous X-ray pulsars. The giant flare events are thought to be due 
■ to magnetic instabilities leading to fra cturing of the crust. The observed QPOs are then associated with shear oscillations 
of the magnetar's crust founcanlll998h . This interpretation ma kes sense since t he fundamental toroidal crustal modes have 
frequencies (~ 30 — 100 Hz) that could match the observations JVan HorrJll98 ph. The properties of crustal modes have been 
invest i gated by several aut h ors, see for example the work by llRm^eiifcCloifir lll98Cll : lMcDermott et all dl98Sh : IStrohmaverl 
Jl99ll) : lMessios et all J200lh : IPirol J2005h . 

In previous studies the basic fact that the magnetic field strongly couples the crust to the fluid core has been ignored. 
This key point, which is wel l known from discussions of pulsar glitch relaxation <|Eassonlll979llAbnev et all99fifr . was recently 
emphasised bv lLevir] J2006h . A rough estimate for the crust-core coupling timescale is provided by the Alfven crossing time, 

t A = 2R/v A ~ 70 B^ 1 p\{ 2 R 6 ms (1) 

where the magnetic field, density and radius are normalised according to B — 10 15 Bi5 G, p — 10 14 pi4 g/cm 3 and R = 10 6 i?6 
cm, respectively. We have also introduc ed the Alfven velocity v\ = B 2 /4np. The cou pling timescale is comparable to the 
period of the fundamental crustal mode jMcDermott et aJll988l Hansen fc Cioffilll980l) . 

P £ ° w 60i? 6 [£(£+ I)]" 172 ms (2) 

This means that, for parameters relevant to the crust-core interface of a magnetar, an efficient coupling to the fluid core 
is already established within a single oscillation of a "crustal" mode. This has crucial implications when one attempts to 
calculate the properties of the modes. In essence, the notion of modes confined to the crust is no longer useful and one is 
forced to consider global oscillations of the coupled crust-core system. 

In this Letter we demonstrate that a magnetically coupled crust-core system admits oscillations with frequencies that 
matches the observational data. Clearly, any initial disturbance in the crust region (say, following a starquake indu ced by 
a mag netic field eruption) will generically shake the field lines and launch Alfven waves into the core. In contrast to iLevinl 
( 200GJ) , we do not think of this process as a damping mechanism. The reason for this is that the Alfven waves in the core will 
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be capable of reaching the crust at the opposite side of the star and be reflected back. By the time that a given wave-packet 
is eventually attenuated by viscosity it would have transversed the core some ~ 10 7 times. Hence, we believe that the correct 
approach to the problem is to consider the coupled crust-core system and compute global oscillation modes. 



2 A SIMPLE TOY MODEL 



We consider a plane-par allel "star" where the fluid core is sandwiched by two slabs of "crust" (this model is similar to the one 
employed bv lPirol 1(2005)). The z coordinate runs from z = +R (the surface) to z — (the core's centre) and ends up back at 
the surface, at z = —R. The crust-core interface is located at z = ±R C . The crust is elastic, with a uniform shear modulus 
p. Furthermore, uniform density, incompressibility and ideal MHD conditions are assumed everywhere. In the unperturbed 
configuration the magnetic field is B = B z. 

The Euler equation in the crustal region is (assuming a harmonic time dependence £ ~ exp(iat)), 



°" ii = - "3-T + 1 — ( V x b) x Bo , 

p dx k 4np 



where 



(Tij = p 



dxi dx 1 



(3) 



Oik and b are the perturbed shear tensor and magnetic field, respectively. Since the fluid is incompressible and all quantities 
vary only with respect to z we have, 

b = B o [*0»e" + yd z t v ] (4) 

Following IPirol J2005T) . we mimic the geometrical factors of the true spherical problem by making the identification 



?+ 1)£,/R . To introduce i in this way is, of course, artificial. In the real spherical problem, the ^-multipoles 



follow from separation of variables. Here they are simply introduced to impose the "correct scaling" for the crust mode 
frequencies. Eqn. J3J now becomes 



[ p, + v\ ] dlii + 



, £(£+1) 

M IP 



0, 



x,y 



(5) 



where p — p/p. Clearly, identical equations describe £ a and £ y , therefore we shall only solve for the former. In addition, as our 
model has reflection symmetry with respect to the "centre" z — 0. Hence it is in p rinciple sufficient to solve only for z > 0. 

T he boundary conditions are formulated in terms of the traction components Jstrohmaverlll99ll : ICarroll et aJll986l IPirol 
2005). For our combined fluid/magnetic field system the traction takes the following form in the crust region: 



T l = p[p + v\]d z C, i = x,y 



and 



T z = 



(6) 



The corresponding expressions for the core region are obtained by setting p — 0. Another important element is the perturbed 
electric field, 

ioBo 



e = ^ x Bo = 



(7) 



c c 

The appropriate MHD boundary conditions are: (i) continuity of the tractions at the crust/core interface, (ii) vanishing 
tractions at the surface and (iii) continuity of the (normal) transverse components of the (magnetic) electric field at the 
interface. These translate into, 



fW) = fTO, d z f(R) = and d,? (Jf) = 1 + P/vl] d z £ x (K) 



(8) 



These constraints imply that the displacement £ is a continuous function at the crust-core interface. Note also that, in solving 
the problem one has two options. One can either solve the full problem and impose conditions analogous to @ for negative 
z. Alternatively, one can use the symmetry of the problem and impose a suitable condition on the solution at the origin (the 
eigenfunction should be either an odd or an even function). 

We begin by considering the special case where the core is decoupled from the crust. As we will see later, this case is 
artificial in the presence of a magnetic field. However, the magnetic coup ling between the crust and the fluid core can b e 
deactivated by altering the boundary conditions at the bottom of the crust JCarroll et ail98d |Pirol2005t iMessios et afeoOll) . 

Solving the Euler equation in the crust provides the displacement, 



£ = Ci e + + C-2 e 



where 



R 



£(£■ 



-21 V2 

a 



1 + v\/p 



and 



R 

V/2 



(9) 



The effective decoupling between the core and the crust can be achieved by requiring vanishing tractions at the interface. We 
stress the fact that this is not the correct interface condition in the magnetic problem. Disregarding this we get, 



Cz 



(10) 
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The requirement of vanishing tractions at the surface gives a similar relation 

C 2 = Ci e 2Ra + 

These two expressions are compatible provided that 



= 1 
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with n = 0, 1, 2, 3, 



(11) 



(12) 



where A = R — R c is the crust thickness. This is th e prediction of our toy model for the "toroida l " crust mo de frequencies. 
It is in good agreement with more rigorous results JlVIcDermott et "al ll988l: lHansen fc Cioflill98d |Piroll2005h . However, the 
model is inaccurate for £ — 1 since a nodeless eigenfunction is allowed, a feature not found in the rigorous spherical problem. 
Given the artificial introduction of £ in the model, this is not surprising. For £ > 2 (and 1 = 1, n > 1) the model predicts the 
correct number of nodes for the eigenfunctions. It is useful to note that the fundamental (n = 0) and first overtone (n = 1) 
quadrupole modes take the respective values cto = 2.45 and o\ = 31.5 (for the latter we have assumed A = 0.1R and va ~ 0). 

In the case of the full crust-core model we solve the Euler equations in the two regions. This (again) leads to a solution 
of the form Q in the crust, while in the core we find 



£ = Die 



i/3. 



+ D 2 e 



03- 



where 



P± = ±<t/va 



(13) 



Once more, the surface boundary conditions enforce (I1H . On the other hand, eqn. (1101 is no longer valid. Instead, the 
appropriate interface conditions result in a 2 x 2 homogeneous system for the coefficients D\ and D 2 . The requirement that 
the relevant determinant must vanish provides the following condition 

n 2 



-2ifl c /5 + 



1 + 



/' 



A 



sinh(Aa+) — i/3+ cosh(Aa+) 



a+ 1 + 



sinh(Aa+) + ij3+ cosh(Aa+) 



which determines the complete spectrum of global oscillations. For the amplitudes we find, 

„2Aa + 



Di 



- Cl lh e 



i + 4 

V A 



1 



-i0+Rc 



and 



Dz 



±Di 



= (14) 



(15) 



(where the upper/lower signs correspond to the same solution). Up to normalisation, this completely specifies the solution. It 
is clear that the modes are either odd or even functions with respect to the origin. It is also worth pointing out that, since the 
limit v\ — > is singular, some care is required in analysing how the pure crust modes emerge as the magnetic field vanishes. 



3 RESULTS: GLOBAL MHD MODES 

Despite being very simple, the toy-model provides a set of interesting results. First note that we can obtain analytic solutions 
to 1141 for a weak magnetic field. When v\ <C p, acceptable mode solutions coincide with the Alfven frequencies: 

oa = (kiv/2R c )vA, k = 0,1,2,3,... (16) 

Remarkably, these frequencies also satisfy 1141 when e 2Aa + = 1, that is, when the frequency coincides with a crustal frequency 
<r„. This triple intersection is naturally interpreted as a resonance between the crust and the core. 

We have determined the exact spectrum by solving l|14|l numerically. A part of the spectrum for £ = 2 and A = 0.1R is 
shown in the left panel of Fig.Q together with the Alfven frequencies l|lri|l . as a function of the ratio 

0.04 B?^ 1 (17) 

We have assumed the value ft ^ 2x 10 16 cm 2 /s 2 jDouchin fc Haenselll200lh . Figure Q shows that, for a given value of the 
ratio 1171 . the spectrum consists of a semi-infinite family of modes (we only show the first few). The separation between 
consecutive modes increases as the ratio attains higher values, i.e. if we increase the magnetic field while keeping the other 
parameters fixed. The most important feature to note in the spectrum is that there are always mode-frequencies comparable 
to the crustal frequencies a n . Fig. Q depicts the spectrum in the vicinity of <to but a similar picture arises for all frequencies 
cr„ of a given I. As is apparent in the figure, for certain discrete values of the ratio v\/jx the mode frequency exactly coincides 
with a crustal frequency a n as well as one of the Alfven frequencies 1161 . For any other realistic value of v\/fl there is always 
a mode with frequency very close (with at most a few percent deviation) to each "crustal" frequency er n . 

The coupled crust-core system obviously has a much richer set of oscillation modes than the decoupled elas tic crust. Yet one 
can in terpret most of the observed QPOs as pure crust modes for various values of I, see the discussion of lstrohmaver fc Wattsl 
( 2005). At first sight our results may seem at variance with this interpretation. After all, our system has a number of additional 
modes. If they are not observed we need to explain why. To do this, we take the standard model for the giant flares at face 
value. The magnetic field erupts and induces a starquake in the crust. On the grounds of physical intuition we would expect 
that the initial perturbation in the crust will predominantly excite those global modes which communicate the least amount 
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Figure 1. Left panel: A part of the global mode spectrum for the toy model of a magnetar with crust-core coupling, see text for details. 
The data corresponds to i = 2 and A = 0.1-R. Note the triple intersection between the exact mode solution from eqn. 1141 . the Alfven 
modes from eqn. 1161 . and the fundamental crustal frequency cto = 2.45. The black dots correspond to observed QPOs frequencies for 
SGR 1806-20, see text for discussion. Right panel: Assessing the excitation of various modes by means of the ratio between the mode 
energy in the crust E crus t and that in the COrG, -C/corc* 

The dashed line labels the crustal frequency cto. 

of energy to the core. To test this idea, we calculate the total energy (kinetic + magnetic) associated with each mode and 
consider the ratio E CIUSt /E COTC . As can be seen from the right panel of Fig.0 the results corroborate our intuition. The mode 
which maximises the energy ratio, i.e. should be easier to excite from an initial crust motion, is the one nearest to the resonant 
frequency a a = a a (analogous behaviour to that depicted in the right panel of Fig. is found for modes in the vicinity of all 
higher a n frequencies, for any I). This would explain why the other modes of the system are more difficult to excite, they are 
predominantly core- Alfven modes which would be energetically more expensive to excite via the proposed mechanism. Similar 
conclusions can be drawn from the displacement £ x , displayed in Fig. [5] for a fixed value of v A /fi, as we move upwards in 
the spectrum. The oscillation amplitude is generally far greater in the core, but the modes located in the vicinity of the crust 
frequencies a n are exceptions. In those cases the crust oscillates with a comparable amplitude. However, these modes cannot 
be considered as "crustal". The eigenfunctions clearly show that they represent global oscillations. 

The conclusion of this discussion is that once the crust is shaken by a starquake the core-crust system will naturally choose 
to vibrate in those global modes which have frequencies similar to the frequencies of the toroidal modes of the uncoupled 
crust, with identical values of £ (and number of nodes in the crust). Hence, the model provides an explanation for all observed 
QPO frequencies, both intermediate (~ 30 — 155 Hz, different values of I) and high frequency (~ 625 Hz, overtone with one 
node in the crust). Our model differs from previous ones only in that the modes are global, not localised to the crust, and 
hence have a significant amplitude in the core. 

A key feature of Figure is the existence of modes with frequencies below the fundamental crustal frequency ctq. This is 
interesting since the presence of low frequency QPOs has been confirmed in the data of both the RXTE and RHESSI satellites 
lllsrael et al|l2005HWatts fc StrohmaveilboOrl) : 18, 26 and 30 Hz QPOs for the December 2004 gi ant flare in SGR 1806-20. It 
is natural to identify the 30 Hz QPO with the fundamental crustal frequency <ro ijlsrael et afeoosh . since the higher frequency 
QPOs then fit the predictions from Eq. J2J quite well. However, we then find ourselves left with a puzzle: what is the origin 
of the remaining two frequencies? There are certainly no toroidal crustal modes with frequency below ao. Our model offers a 
natural explanation for these low-frequency modes. 

Although our model is not designed to provide us with quantitatively accurate predictions, we can still attempt to match 
the observed data for the low frequency QPOs of SGR 1806-20. To do this we first identify cr = 2.45 with the 30 Hz 
frequency (presumably, this identification can be made rigorous in a calculation for a realistic neutron star model). Then the 
frequencies 18 and 26 Hz correspond to modes with a — 1.46 and 2.12, respectively. Some navigation in Fig. leads us to 
the value v\ — 0.055/i which has all three desired modes (indicated by filled circles in the figure). This value is reasonably 
consistent with Bi$ ~ 1 assuming P14 sa 1 ( note that the magnetic field estimated from the spin-down rate for SGR 1806-20 
is 7.8 x 10 14 G \ Woods fc Thompson! 120041) . Moreover, the energy argument indicates that it is reasonable to expect these 
QPOs to be excited, cf. Fig.Q Finally, we have an additional mode at « 22 Hz. In a sense, this could be considered a testable 
prediction. It would certainly be rewarding if a QPO with this intermediate frequency were to be found in the data! 

4 CONCLUDING REMARKS 

We have used a simple plane-parallel toy-model to discuss global MHD modes of a neutron star. The model, which takes 
into account the magnetic coupling between the elastic crust and the fluid core, provides results that are highly suggestive. 
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Figure 2. Fluid displacement £ x for £ = 2, A = 0.1R and v\ = 0.055 jl. The modes illustrated are, from left to right: a = 0.74, 2.45, 9.87. 
In the middle panel note the comparable oscillation amplitudes in the crust and the core for the mode with a m = a\. 

The system is characterised by a rich spectrum of global modes "living" both in the crust and core, with a mixed magneto- 
elastic identity. Among myriads of available modes (for a given multipole £ and ratio of Alfven velocity to shear velocity) the 
system will naturally excite the mode(s) which transfer as little energy as possible to the core, thus minimising the overall 
energy budget. We demonstrated that these favoured global oscillation modes have frequencies that are similar to the "pure 
crustal" frequencies from the non-magnetic problem. Our model provides support for the asteroseismology interpretation for 
the observed magnetar QPOs. In addition, it provides a natural explanation for the presence of lower frequency (< 30 Hz) 
QPOs seen in the December 2004 giant flare of SGR 1806-20. 

Even though our model is based on drastic simplifications, the basic physics predicted should survive for more realistic 
neutron star models. Obviously, one must expect corrections (and extra technical complications!) once the real problem is 
considered. A more complex magnetic field configuration will lead to intricate coupling between different multipoles, non- 
uniform density/crust elasticity and core stratification leading to buoyancy effects are likely to affect the mode structure, 
and the expected superconductivity will affect that nature of the Alfven waves. Nevertheless, global MHD modes due to the 
inescapable crust-core coupling are still likely to emerge, leading to intricate resonance phenomena. 
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